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Abstract: 

In this study, experiments were performed to 
determine parameters effecting Material Removal rate 
(MRR) along with performance measurement analysis 
with respect to respect to input machining parameters. 
Experimental work was conducted on EN-19 alloy steel 
with Copper as tool electrode and EDM oil as dielectric 
fluid. The data observed during experimentation has 
been used to yield responses in respect of material 
removal rate (MRR). The objective of this paper is to 
study the influence of operating input parameters of 
copper electrode on material removal rate of EN-19 
material followed by optimization for confirmation test 
purpose. The effectiveness of EDM process with 
copper electrode is evaluated in terms of the material 
removal rate. In this work the influence of the 
parameters such peak current, voltage gap, pulse on 
time, duty cycle and flushing pressure. Material 
removal rate (MRR) in this experiment was calculated 
by using mathematical method. The result of the 
experiment then was collected and analyzed using 
Minitab software. This was done by using the response 
surface method and Anova analysis. 

Keywords : Material Removal Rate, Response Surface 
Methodology, Regression. 

I. INTRODUCTION 

Electric discharge machining has extensive 
applications for manufacturing dies and tools to 
produce mouldings, die casting, and sheet metal dies 
etc[1][2]. Implementation of EDM process will awaken 
manufacturing engineers, product designers, tool 
engineer and metallurgical engineers about unique 
capabilities and benefits of this process[3].EDM can 
be used for machining of high precision of all type of 
conductive material (metals, alloys, graphite, 
ceramics etc.) of any hardness. In EDM process [4], 
material removal from work piece is done by means 
of a series of electrical discharges. This paper 
presents work on machining by EDM for EN-19 alloy 



steel. Certain parameters in EDM process directly 
influence the process outputs. Setting appropriate 
values for such parameters requires the 
implementation of many machining trials. This leads 
to time consuming and expensive experimental work. 
Response Surface Methodology (RSM) has been 
used for modelling EDM machining of rectangular slot 
size 15 mm x 20 mm on EN-19 material using copper 
electrode tool [5, 6-9]. Response surface method is 
employed to signify relationships between inputs and 
significant outputs based on minimum number of 
experiments. This paper presents a mathematical 
modelling of EDM machining process on EN-19 alloy 
steel using RSM approach. 

II. EXPERIMENTAL WORK 

The material used for this work is EN-19 alloy steel 
square plate of size 100mm x 100mm x 20mm with 
density 7.85 g/cm 3 . The specimen is machined on 
conventional milling at depth of cut 0.25 mm to 
produce a plane surface. Copper electrodes (99.97% 
pure, density 8.96 g/cm 3 and melting point of 1086 
•C), parallelepiped shaped 20mm x 15mm 100mm is 
used in the experiment. The machine used is ENC 
EDM Microcut make with NC control in Z-direction 
with EDM oil as dielectric medium. 

III. EXPERIMENTAL PLAN 

Experiments are planned on the basis of RSM 
technique used in experimental design. The codes 
are calculated as functions of the range of interest of 
each factor as shown in Table 1. A central composite 
design with five input variables having five levels 
between (• = ±2) coded values and 32 experimental 
runs were performed. Different variables represented 
by x v x 2 x 3 x 4 x s and their levels are given in Table 
2. The coded numbers for the variables used in tables 
are obtained from the following relationship[10] : 
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Table 1. Coded & actual values[10] 

The numbers of test required are chosen with the 
standard 2 k full factorial central composite design. 
CCRd provides as much as information as a five level 
factorial, requires many fewer tests and has been 
shown to be sufficient to describe the majority of 
process responses [11, 12, 13]. Each experiment is 
performed using copper electrode, with a particular 
set of input parameters chosen randomly from the 
planned set of experiments. The polarity of the 
electrode is set as negative. The depth of machining 
is set at 2mm for all sets of experiments. 
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Table 2. Coded levels 

IV. REGRESSION MODELLING & ANALYSIS 

According to the experimental plan a total of 32 
experiments are conducted, each having the 
combination of various values of process variables x,, 
x 2 x 3 x 4 x s Each of the responses is fitted into a linear 
equation represented by: 

Y (MRR) = -0 +;x, + > 2 x l + - 3 x 3 + - t x 4 + - s x^~ (1) 

Where, Y is the response and x v x 2 x 3 x 4 x s are 
coded levels of the variables. The coefficients -0, ■, - 2 
• 3 •„ - 5 can be calculated by solving the following 
equation: 



= (x T x)VY- 



(2) 



where, • is the matrix of parameter estimates, x is 
the matrix of independent variables, x T is the 
transpose of X matrix and Y is the matrix of measured 
responses. Table 3 gives the design matrix and the 
responses. Analysis of variance (ANOVA) is 
performed to test the adequacy of the proposed 
models. The variance ratio denoted by F in ANOVA 
tables, is the ratio of the mean square due to a factor 
and the error means square. In this robust design F 



ratio can be used for qualitative understanding of the 
relative factor effects. A large value of F means that 
the effect of that factor is large compared to the error 
variance. So the larger value of F, the more important 
that factor is in influencing the response [14]. 

In this work from Table 4, Anova table shows the 
most important factor is input current with 311.69 F 
ratio and pulse on time with F = 6.59. F ratio of other 
factors is not significant and has minimum effect. 
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Table 3. Design matrix and Response 

Anova is used to test the null hypothesis with 
regard to the data gained through experiments 
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[16]. Empirical models are fitted for the stated rate is given by Eq. (3) and its analysis of variance is 

responses material removal rate. Analysis of variance in Table 4. 
is carried out on all the fitted models for a confidence 

level of 95%. The fitted model of material removal MRR = 6.8 + 0.202 T + 0.401 V mn + 1.48 DC + 



2.79 I + 0.2 F p (3) 



Source 
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Table 4 Anova of MRR 



From Equation 3, the factors gap voltage, input 
current, duty cycle, pulse on time have an additive 
effect on the material removal rate where as flushing 
pressure has minimum impact on MRR. Analysis of 
the residuals of the model shown in Equation 3 is 
performed to test assumptions of normality, 
independence and constant variance figure 1 of 
residuals. The quantitative test methods mentioned 
earlier are employed again, and none of the 
assumptions are violated. 



Residual Plots for MRR 
Normal Probability Plot Versus fits 




Figure 1. Residual plots for MRR 

Regression analysis is carried out to ensure a least 
squared fitting to error surface in Minitab 15 
environment. Regression analysis has been 
performed to find out the relationship between input 
factors and MRR. During regression analysis it is 
assumed that the factors and the response are 
linearly related to each other. The general first order 
model is proposed to predict the surface roughness 
over the experimental region can be expressed as 
Equation 1 . 

In general, the R 2 adj statistic will not always 
increase as variables are added to the model. In fact, 
if unnecessary terms are added, the value of R 2 adj 
will often decrease. When R 2 and R 2 adj differ 



dramatically, there is a good chance that no 
significant terms have been included in the model 
[15]. For this experiment the R 2 value indicates that 
the predictors explain 95.50 % of the response 
variation. Adjusted R 2 for the number of predictors in 
the model was 91 .0 % both values shows that the 
data are fitted well. The prediction model was then 
validated with another set of data. The verification of 
the tests results for material removal rate is shown in 
Table 5. The predicted machining parameters 
performance is compared with the actual machining 
performance and a good agreement is observed 
between these performances. In Table 5 process 
factors are given in terms of natural factors and their 
corresponding coded factors. In order to assess the 
accuracy of the prediction model, percentage error 
and average percentage error were recorded. 
Percentage of prediction errors is shown in the last 
column of Table 5. The maximum prediction error was 
17.38 % and the average percentage error of this 
method validation was about 6.87%. As a result, the 
prediction accuracy of the model appeared 
satisfactory. 
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Table 5. Error Prediction 
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Surface Plot of MRR vs Ton, Ip 




Figure 2. Surface Plot for MRR Vs Ton & Ip 



Surface Plot of MRR vs Vgap, Ton 
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Figure 3. Surface Plot for MRR Vs Ton & Vgap 



Surface Plot of MRR vs Vgap, Ip 




Figure 4. Surface Plot for MRR Vs Vgap & Ip 

Figure 2 show that material removal rate increases 
with the increase in input current. Though, there is 
low MRR at the low I and low T on which gives direct 
linear relation but there after it is increased with rise in 
current and pulse on time. Figure 3 show similar 
relationship between MRR versus T on and V gap in 
which it clear that there is linear relation of V„ a „ with 

9 a P 

MRR. The influence of F p is minimum on material 
removal rate. And Figure 4 show the relation between 
l p and Vgap. In this graph the high peaks show that 
the there is increase in MRR with increase in V„ a „ 



V. OPTIMIZATION OF MRR 

The designed experiments involve determination of 
optimal conditions that will produce the "best" or 
"optimum" value for the response (MRR). Depending 
on the design type (factorial, response surface, or 
mixture), the controllable operating conditions may 
include one or more of the following design variables: 
factors, components, process variables, or amount 
variables. Optimal settings of the design variables for 
one response may be far from optimal or even 
physically impossible for another response. Response 
optimization is a method that allows for compromise 
among the various responses. The optimization is 
carried by obtaining the individual desirability (d) 
for each response combining the individual 
desirability to obtain the combined or composite 
desirability (D) thereby maximizing or minimizing 
the composite desirability and identifying the optimal 
input variable settings. Here in case of surface 
roughness optimization, it single response 
optimization where the overall desirability is equal to 
the individual desirability. 

V (a) INDIVIDUAL DESIRABILITY 

As in this case of MRR we need to optimize single 
response, so here individual desirability (d) for 
material removal rate is obtained using the goals and 
boundaries for MRR that is given in Minitab session 
window. There are three optimization goals desired 
as follows: 

• minimize the response 

• target the response 

• maximize the response 

For material removal rate (MRR) it is desirable to 
obtain maximum value for better surface finish of 
material. As response MRR is desired to be 
maximizing for which determination of target value 
and an allowable maximum response value is 
provided to response optimizer. The desirability (d=1) 
is one for MRR response below the target value: 
above the maximum acceptable value the desirability 
(d=0) is zero. 

In the below Table 6, y is the response value, T and 
L are the target and lower boundaries (i.e. minimum 
and maximum acceptable values for the response), 
respectively, and T is the target. For the MRR(y) to 
maximize by: 
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f,(y) = 
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Y< L 


y-L/ T-L 


L • y • T 


1 


Y > T 



Table 6 Maximization of response by individual desirability 



5(b) RESPONSE OPTIMIZATION 



Lower Target Upper Weight Import 



| MRR Maximum 85 175 175 0.1 1 



Table 7 EWR Range 

Starting Point: Ton = 50; Vgap = 12; DC = 5; Ip = 25; 
Fp = 0.3 



Ton = 


90 


Vgap = 


18 


DC 
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IP = 


45 


Fp = 


0.5 



Table 8. Global Solution 

Predicted Responses as shown in Minitab 
session window. 

MRR = 168.053, desirability = 0.992000 
Composite Desirability = 0.992000 
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Figure 5. Optimization Plot 

Each response in the research work are expressed 
separately as linear and non linear functions of input 
variables such as I T on , V gap , DC, F . Now aim is to 
maximize the response MRR and simultaneously 
maintain other responses in EDM process. As shown 
in Table 8 global solution of input parameters is 
obtained by response optimizer. To determine global 



solution of input variables in order to satisfy the above 
criteria of MRR maximization, it had been solved by 
Response optimizer desirability maximization function 
in Minitab 15 environment. The individual desirability 
for MRR material removal rate is 1. To obtain this 
desirability, the optimum values factor levels can be 
set as shown under Global Solution in the Minitab 
Session window in Table 8. That is, \ = 45, Ton=90, 
V gap =18, DC=7, F p = 0.5. The optimum predicted value 
for MRR = 168.053 obtained for 99.2 % desirability. If 
it is further desired to improve this initial solution, you 
can use the plot. Move the red vertical bars to change 
the factor settings and see how the individual 
desirability of the responses and the composite 
desirability change. 

VI. RESULT AND DISCUSSION 

Results show that input current, pulse on time, duty 
cycle, voltage gap are significant factors for MRR and 
flushing pressure has minimum effect on the material 
removal rate of EN-19 alloy steel material. Finally, a 
mathematical model was developed using multiple 
regression method to formulate the input current, gap 
voltage, and pulse on time, and flushing pressure to 
the MRR. 

VII. CONCLUSIONS 

From this analysis we can conclude that the most 
significant EDM process variable influencing all the 
stated machinability parameters of EN-19 alloy steel 
is pulse current. The significance order of other 
parameters is gap voltage followed by pulse on time 
and gap voltage. These models can be effectively 
utilized by the process planners to select the level of 
parameters to meet any specific EDM machining 
requirement EN-19 alloy steel within the range of 
experimentation. 
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